How much of one-way computation is just thermodynamics? 
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In this paper we argue that one-way quantum computation can be seen as a form of phase 
transition with the available information about the solution of the computation being the order 
parameter. We draw a number of striking analogies between standard thermodynamical quantities 
such as energy, temperature, work, and corresponding computational quantities such as the amount 
of entanglement, time, potential capacity for computation, respectively. Aside from being intuitively 
pleasing, this picture allows us to make novel conjectures, such as an estimate of the necessary 
critical time to finish a computation and a proposal of suitable architectures for universal one-way 
computation in ID. 



I. MOTIVATION 

In a thought-provoking paper [l[, Toffoli argued that 
physics is computation. He achieves this by mapping 
the Lagrangian 'action' in physics, which counts in how 
many different ways a system can evolve dynamically, to 
Shannon's entropy in information theory, which quanti- 
fies the lack of information one has about the state of a 
system. So, "action is to dynamics exactly what entropy 
to statics" [l|. In this paper we will invert Toffoli's logic 
and map the process of a one-way quantum computa- 
tion 01 to a physical description in the language of ther- 
modynamics. Ours is a conceptual and heuristic paper 
whose main aim is to put one-way computation in a much 
broader physical framework. An exciting consequence of 
this idea will be that one-way quantum computation can 
be viewed as a phase transition from a disordered phase, 
with diverse possible outcomes of the computation, to an 
ordered phase, where the unique solution to the compu- 
tation is reached. 

Adopting this view has a number of merits. First of 
all, some already known results in one-way quantum com- 
putation become transparent when formulated with the 
view that the number of ways a computation can be im- 
plemented determines computational capacity. For ex- 
ample, the fact that ID one-way cluster computer is not 
universal as well as the fact that the 2D cluster computer 
can be. Secondly, it allows us to explain some desirable 
properties of the one-way computer design without hav- 
ing to use any heavy mathematical machinery. Here we 
will address the issue of the efficiency of one-way com- 
puters as a function of their inter-connectivity. 

One-way quantum computation relies on the initial, 
highly entangled state of many qubits, for instance the 
cluster state, which forms the resource for the computa- 
tion. An algorithm can be implemented by measuring a 
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fraction of the state's qubits in a specific and controlled 
way thereby driving the rest of the 'computer', i.e. the 
state of all qubits, to produce the solution to the algo- 
rithm. Since the computation is driven by measurements 
and not by unitaries, the cluster state computer is not 
reversible, which is why it is referred to as one-way com- 
puter [TtI . The one-way computer can only work by way 
of the entanglement that correlates qubits inside the clus- 
ter state and 'spreads' information about the performed 
measurements over all qubits in the computer. Entangle- 
ment is thus the true resource of the computation which 
is used up to produce the 'solution', i.e. bits of informa- 
tion. It is, in this respect, analogous to solid state sys- 
tems that trade energy for reduced entropy when they 
undergo a phase transition. 

Moreover, the one-way model of computation is fully 
equivalent in power and capacity to the standard model 
of a quantum Turing machine Understanding the 
key ingredients and the dynamics of one-way computing 
by formulating it as a critical physical process will shed 
light on what the power of quantum computing is, that 
goes beyond the realm of classical computing. Further- 
more, the measurement-based aspect of one-way comput- 
ing offers practical advantages over the conventional gate 
model in a number of experimental settings making it the 
model of choice for quantum computing. A recent exper- 
iment has already demonstrated the first implementation 
of basic gates and the Grover search algorithm using en- 
tangled photons as cluster states [1]. Finally, discussing 
the one-way quantum computer as a thermal system will 
help identifying new candidates for realistic experimental 
implementations . 

The paper is structured as follows. We summarise the 
basics of one-way computing with the cluster state as 
the resource in Section [TTl In Section IIIII we briefly ex- 
plain the idea behind the concept of phase transitions 
and review Peierls' argument of why there is no phase 
transition in the ID Ising model and why there can be in 
higher dimension. We present our main ideas in Section 
IIVI where we propose a mapping of relevant thermody- 
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namical quantities to central computational characteris- 
tics, in particular the free energy to the computational 
potential. In Section |V] we derive first results from the 
established heuristic analogy and give an outlook on fur- 
ther questions that could be addressed using the same 
link. 



II. ONE-WAY CLUSTER COMPUTING 



state, which has the mathematical form 

Ibox-clustcr) = |0)|+)|0)|+) + |0)h)|l)h) 

+ |l)h)|0)|-) + |l)|+)|l)|+). (3) 

Pictorially the box cluster can be represented as a square 
with one of the four qubits at each vertex. 



A possible initial state of the one-way computer, the 
cluster state, is prepared in the following way. First take 
an array of N qubits arranged in a d-dimensional square 
lattice. Then prepare each of these qubits in an equal 
superposition state with respect to the z-basis, |0) and 
as \+) = ^(|0) 4- Finally, apply a two-qubit 

C- phase gate, CZ, to every pair of nearest neighbour 
qubits in the lattice, changing, for instance, | + +) into 
the entangled state 
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The full two-qubit C-phase operator can be written in 
matrix form with respect to the basis |00), |01), |10), |11) 

as 
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The resulting fully C-phased, highly entangled A^-qubit 
state is the cluster state which can be depicted schemat- 
ically as below for A = 16 in the two-dimensional case. 




The red dots represent the qubits and the lines connect- 
ing them signify that they have been C-phased. Note, 
that the order of applying the C-phase gates does not 
matter since the C-phase operators commute. Experi- 
mentally such an operation could be achieved by letting 
an Ising-type Hamiltonian act on the whole lattice for an 
appropriately chosen time. 

The amount of entanglement in the cluster state is di- 
rectly proportional to the number of applied C-phase 
gates [l^. A suitable measure of entanglement is the 
relative entropy of entanglement The relative en- 

tropy of entanglement of a state p is defined as Efj{p) = 
mino-gu ^(/allfT), where S'(/ci|j(T) = tr[/9 log p—p log cr] is the 
quantum relative entropy and the minimisation is taken 
over all disentangled mixed states a. To illustrate that 
the relative entropy of entanglement is a useful measure 
in the context of cluster states consider the box cluster 



All cluster states of any size are pure states and their 
relative entropy reduces to S{p\\(t) = — tr[p log a]. The 
relative entropy of entanglement then requires to find 
the disentangled mixed state a that minimises the rel- 
ative entropy. It can be shown that this minimum 
is reached for a = diag{ , . . . , } so that the rela- 
tive entropy of entanglement of the box cluster becomes 
I box-cluster)) = 4. This coincides with the number 
of applied C-phase gates. The same holds for general 
cluster states where the entanglement is always propor- 
tional to the number of applied C-phase gates Q . 

To see the mechanism of one-way computing at work, 
consider the simplest case of a two qubit cluster state 
computer (see [2| as the general reference). Instead of 
being prepared in |-|-), the first qubit is the input and 
can be any state \tp) = a\0) + /3|1), with a and (3 arbi- 
trary complex numbers with |Q;p-|-|/3p = l. The second 
qubit is initially prepared in the state |+) and will serve 
as the output state of the computation. The two qubits 
are C-phased and become entangled in the total state of 
both qubits, q;|0)|+) +/3|1)|— ). As a computational step 

we now apply a Hadamard gate, i7=--^^| |^,to 

the first qubit which leads to new total state q;|+)|+) + 
/?|— )|— ). Rewriting this state we find that it can be cast 
in the form ■^[\0)(g)H\iJj) + \l)'»a:, H \ip)], with being 
the usual Pauli matrix. This implies that a measurement 
of the first qubit in the computational z-basis will in any 
event have the effect of tclcporting the input state to 
the second qubit - up to a Pauli transformation. Indeed, 
the state of the second qubit after measuring the first 
qubit can be written more compactly as cr™ H \ tp) where 
we have introduced the new variable m £ {0, 1}. The 
variable m corresponds to the two possible outcomes of 
our measurement in the computational basis. After mea- 
suring the first qubit and finding either one of the two 
values, we have to correct and apply the inverse Ha™ to 
obtain the fully teleported state lip). 

This simple example can be generalised to implement 
any rotation Rz{(j>) of the arbitrary input qubit where 
the angle of the rotation is varied by using the basis 
{--i=(|0) ± e"^|l))} for measurement. This will leave the 

second qubit in the desired rotated state, al^HR^icf)), 
again up to a Pauli transformation. Furthermore, any 
single-qubit rotation U £ SU{2) can be decomposed as 
U = Rz{l) RxiP) Rz{a), where a, j3 and 7 are the three 
Euler angles of the arbitrary rotation. To be able to im- 
plement any single-qubit unitary using one-way cluster 
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computing we require a linear 4-qubit cluster. With the 
same method as before the final unitary implemented on 
qubit # 4 is 

U = <i Hi?,((-l)™^7)i?,((-l)"i/3)i?,(a). 

(4) 

This form reveals most clearly how the subsequent mea- 
surements in the one-way cluster computation have to 
be adapted. The angles of the rotations and therefore 
the final correction depend on the outcomes of previ- 
ous measurements, specified by the variables mi, TO2 and 
7713, which imposes a time ordering in the whole scheme. 
Equivalently to applying active single-qubit operations 
such as rotations and then measuring in the computa- 
tional z-basis, one can immediately include the single- 
qubit operations into the measurement process and al- 
low instead for arbitrary single-qubit measurements in 
the X — y plane. 

The above method lies at the heart of one-way compu- 
tation and shows that local, single-qubit measurements 
arc sufficient to manipulate the whole cluster. More ex- 
plicitly, any computation in the one-way computer is 
driven by computational measurements like the above, 
i.e. arbitrary single-qubit measurements in the x — y 
plane specified by the angle cj). The binary outcome of the 
computational measurements is always random and the 
result is fed forward to determine, together with the ac- 
tual algorithm of the implemented computational prob- 
lem, which qubit to measure next and in which mea- 
surement basis (p. All single-qubit measurements that 
can be implemented simultaneously form an 'equal-time' 
set Qtjl^. In general there are a number of such sets 
that needs to be applied in a particular successive or- 
der thereby creating a time-arrow in the computational 
process and we may assign one time-unit to each equal- 
time measurement set Qt- Manipulating the cluster state 
by measuring parts of it will affect the rest of the cluster 
through the entanglement that correlates the parts. This 
process is what makes the computer compute, see lower 
schematic in Fig. [31 

After a sequence of measurements, at time tend, the 
remaining, so far unmeasured qubits of the cluster state 
are ready for readout. Readout measurements arc simple 
single-qubit projections in the z-basis and can be imple- 
mented at one instance in time altogether. The results 
of the measurements after tend are fully deterministic, 
i.e. for any given input, the readout measurements will 
produce a particular bit-value without ambiguity. With 
every deterministic readout measurement we thus find 
one bit of information and the set of all of them presents 
the correct solution of the initial computational prob- 
lem with total information /max- However, if we were 
to stop the computational phase prematurely, this would 
result in a probabilistic scenario where only part of the 
final output may already be deterministically available, 
i.e. some qubits have become disentangled from all other 
qubits and are either in state |0) or in |1) with respect to 
the z-basis, while other qubits are still in superpositions 
with respect to the z-basis thereby making the results of 



their readout measurements probabilistic. In the graphic 
on the right in Fig. [5] we depict the rising information 
I{t) (i.e. deterministic knowledge of bits of the outcome) 
with proceeding time assuming that we have already per- 
formed some equal-time measurement sets Qt up to time 
t and then read out. There are in fact two critical times. 
The first critical value tcrit appears, when our informa- 
tion about the final result becomes non-zero, i.e. when 
we start to gain partial knowledge of the outcome-bits by 
reading out. However, we have to wait until tend where 
the full algorithm has been implemented to obtain the 
complete output, /max- 

We will proceed to compare one-way cluster compu- 
tation to a thermodynamic phase transition process, the 
occurrence of magnetisation in the Ising model. Before 
we go into the analogy let us here first summarise the 
concept of a phase transition and restate Peierls' argu- 
ment why there can be none in the ID Ising model. 



III. PHASE TRANSITIONS AND PEIERLS' 
ARGUMENT 

A phase transition Q is a macroscopically observable 
effect which occurs when a parameter of the system is 
tuned across a critical value, beyond which the system 
takes on a qualitatively different phase. For instance, 
an ordered phase could be a magnetised piece of metal 
or Bose-Einstein condensate (BEC) in a Bose-gas. The 
notion of phase transition is intimately related to that 
of criticality which is a concept of central importance 
in solid state physics. Macroscopic objects, in the form 
of solids, liquids and gases, undergo a diverse range of 
phase transitions under variation of the temperature, or 
external field, or pressure for instance. 

When the two-dimensional Ising spin lattice under- 
goes a phase transition at its critical temperature, this 
means that above this temperature, the spins were in 
a randomised, disordered state, where all directions are 
equally likely, but below that temperature all the spins 
align and point in some (randomly chosen) direction, see 
upper schematic in Fig. [3l This transition from a dis- 
ordered to an ordered state is signified by spontaneous 
symmetry breaking of the z-direction invariance of the 
state of the Ising lattice. An additional parameter has 
to be introduced to specify which particular 'order' was 
chosen by the system. The order parameter in this case is 
the total magnetisation, which is finite below the point of 
criticality and becomes suddenly zero at and above this 
point. Symmetry breaking is, therefore, another indica- 
tor of phase transitions. 

The main advantage of the concept of criticality is the 
fact that wide ranging systems behave very similarly to 
each other at the critical point. This realisation leads 
to the concept of universality, namely that correlations, 
and consequently other observable system properties can, 
when critical, be described using only a small number of 
parameters, that are completely independent of the na- 
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ture of the system. As a result, one can derive some 
rather general conclusions about the existence of phase 
transitions. It is well-known that there are no discrete or 
continuous order parameter phase transitions in ID, no 
matter what particular system we are talking about, pro- 
viding that the interactions are short ranged. This result 
has been proven using many different methods, first by 
Peierls (and tidied-up by Griffiths) Q, then by Mermin 
and Wagner Q , and finally by Hohenberg [l3| • The gist 
of all these arguments is that any thermal fluctuations 
in ID are enough to destroy order that is generated by 
short range interactions. In other words, entropy always 
dominates energy in ID. In this paper we will rely on the 
argument due to Peierls since it is the most suitable one 
for our application to one-way computation. 

The intuitive picture of Peierls' argument [H] is as fol- 
lows. The classical one-dimensional spin chain interact- 
ing with the Ising short range coupling would undergo a 
phase transition from a disordered phase to an ordered 
phase if a non-zero magnetisation arose below a certain 
critical temperature. To test whether the ordered state 
can be an equilibrium state of the system for any finite 
temperature, we assume that the system is initially or- 
dered and then perturb the state slightly while observing 
the change of the free energy, AF, with 



F = U-TS, 



(5) 



where U is the internal energy, T the temperature and 
S the entropy of the state. If the ordered state was an 
equilibrium state, then its free energy must be minimal 
and any perturbation out of this state would increase 
free energy, AF > 0. If, on the other hand, there ex- 
ists a perturbed configuration that can decrease the free 
energy, the system will tend towards this new configu- 
ration preferring it energetically over the ordered state. 
This immediately implies that any order will be destroyed 
and no criticality is exhibited. 

The Hamiltonian of the ID Ising spin chain with near- 
est neighbour interaction is given by 



N 



(6) 



where the aj are the Pauli spin operators in z-direction 
for the j-th qubit and N is the number of qubits in the 
chain. The interaction J is positive in the ferromagnetic 
case and negative in the anti-ferromagnetic case. The or- 
dered state is reached when all spins are aligned, either 
in the positive or negative z-direction. We proceed by 
applying a single random flip, which realises the smallest 
possible perturbation. The change of internal energy un- 
der such perturbation is AU = (i?)flipped — (-ff)ordorcd = 
—J{N — 2) — {—JN) = 2J, because every spin is con- 
nected to two nearest neighbours. There are N different 
ways to flip a single spin so that the change in entropy is 
AS* = kB In TV. Therefore the change in free energy is 
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FIG. 1: How to flip spins in a 2D lattice in the best case 
scenario. This figure shows a possible flipping pattern that 
cuts the lattice into domains of equal alignment. A path 
that wiggles back and forth inverts a lot of nearest neighbour 
couplings and the energy of the lattice increases considerably 
(maximally). However there are many similar, albeit differ- 
ent such paths so that the entropy of the lattice increases at 
the same rate as the energy. Therefore, there exists a finite 
temperature range, where the entropic fluctuations are more 
costly than the energy gain and the fully aligned (ordered) 
configuration remains the minimum of free energy. Above 
this critical temperature, the ordered phase is destroyed and 
a phase transition occurs. 



and since we are looking at the thermodynamical limit 
of large TV it is clear that AF will always be negative 
for any finite temperature. The conclusion is thus, that 
in ID no ordered state can exist in a finite temperature 
equilibrium. The entropy, i.e. the thermal fluctuations, 
will always outweigh the tendency of energy to create 
order. Similar conclusions can be drawn for all short 
range interactions completing Peierls' proof that no phase 
transition takes place in one dimension. 

This situation is changed dramatically in two and 
higher dimensions. In two dimensions the ordered state 
is also the one where all the spins point in the same di- 
rection. However, flipping a single spin is now no longer 
enough to destroy this order; instead the smallest pertur- 
bation must cut the lattice into two different domains [20|. 
Let us now calculate the best case scenario for such a cut 
to destroy order. As far as the balance between change in 
energy and entropy is concerned it is the one which cuts 
the lattice into two halves by crossing as many nearest 
neighbour couplings as possible along the way. For such 
a cut, the change in energy is at most 2A''J, since this is 
the number of couplings present in the lattice (see Fig. 
[T]for an illustration). The entropy change, on the other 
hand, is roughly fc^ In 3^ because after each spin we have 
a choice to proceed with the cut in at most three possi- 
ble directions (everywhere but backwards). The change 
in free energy for this kind of perturbation is thus 



AF = 2N. J - keTN \n3. 



(8) 



By setting AF = we obtain the following critical tem- 
perature for the stability of the ordered phase 



AF = 2J-kBT\nN, 



(7) 



2J 



ks In 3 



(9) 
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This is astonishingly close to the exact temperature ob- 
tained by Onsager [111] (T = 2.27^) especially given 
that the above argument is very simplistic. 

To summarise, in two dimensions the two competing 
quantities can be balanced and a phase transition can 
occur at a finite temperature. With the same argument 
one can show that this holds also true for higher dimen- 
sions. We will now proceed to draw the analogy between 
the phase transition in the Ising model and the process 
of one-way cluster computing. 



IV. THE COMPUTATIONAL ANALOGY 

A central question for one-way quantum computation 
is, 'Why does it work?' Although it is possible to go 
through each mathematical step, using the rules of one- 
way computing, it is not obvious how the computational 
power arises. Clearly, the existing correlations in the 
cluster state have to be traded somehow for the informa- 
tion about the solution of the computational problem. 
But how can we find the optimal way to do it? And 
what is the ingredient that makes, apparently, one-way 
quantum computers more powerful than their classical 
counterparts? 

In the following we compare the process of one-way 
computation to the phase transition in the Ising model 
and identify potentially useful quantities for the one-way 
computer in analogy to those used in thermodynamics, 
as shown in the table in Fig. [51 Our comparison can- 
not directly answer any detailed questions yet. However, 
the ideas sketched here open up a completely new way of 
looking at the dynamics of the one-way computer. Us- 
ing the language of phase transitions it becomes intuitive 
why one-way computing can work and it also sheds light 
on how the computational process can be traced with 
only a few characteristic variables. The analogy will fur- 
thermore allow us to give rough estimates of minimal re- 
quirements for universal quantum computing. One result 
we can re-produce is to heuristically understand why the 
ID cluster state is not a universal resource for one-way 
computation whereas the 2D cluster state has this po- 
tential. Additionally, we will present an estimate of the 
critical time for a computation depending on the dimen- 
sion or connectivity of the resource state for dimensions 
higher then ID. 

We saw how cooling a magnetic material of a multi- 
tude of spins means that the individual spins are flipped 
to obtain a new thermal equilibrium between two quanti- 
ties, the energy and the entropy. Thermal equilibrium is 
determined by the second law of thermodynamics as the 
configuration that maximises the entropy for fixed inter- 
nal energy or, to put it differently, the configuration that 
achieves the lowest possible free energy for a given tem- 
perature. We propose that an optimal one-way computer 
should work in principle in the same way as nature. Sim- 
ilarly to thermal systems, we expect a principle of least 
action that balances two quantities at each step in the 
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One-way computing 






Energy U 
Entropy S 
Temperature T 
Free Energy F 
Cooling 

Magnetisation AI 


Entanglement E 
Computational Capacity C 
Inverse Time j 
Computational Potential P 
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Extracted Information / 



FIG. 2: Tliis table compares the main quantities relevant in a 
thermodynamical phase transition process and in the one-way 
computer model. For details please refer to the text. 



computation to maintain enough computational power 
to finalise the computation. These two quantities are the 
entanglement and the computational capacity which will 
be defined in the following. 

While in the Ising model each spin interacts via Heisen- 
bcrg coupling with its nearest neighbour, in the one-way 
cluster computer each pair of nearest neighbouring qubits 
is C-phased, see Fig. [31 Even more, this C-phase gate op- 
eration can be achieved by letting precisely the Heisen- 
berg interaction act on the qubits for a certain amount 
of time. In the Ising model, the Hcisenbcrg interaction 
results in a non-trivial internal energy U. In thermody- 
namics the internal energy is the resource for a physi- 
cal system to perform work and, in the context of phase 
transitions, energy provides the capability to create or- 
der. It is appealing and intuitive to view entanglement 
in much the same way for the one-way computer. The 
entanglement, which is generated by the C-phasing, 
here allows the transport of information across the cluster 
state, eventually resulting in the ordering of the possible 
outcomes of a computational task into a single, unique 
solution. It is thus the resource for executing dynam- 
ics in the one-way computer. Each spin flip in the Ising 
model costs one unit of energy, similarly each computa- 
tional measurement in the one-way computer consumes 
approximately one unit of entanglement. The number of 
gates implemented in this way, i.e. the computational 
'work', is proportional to the number of units of entan- 
glement used up by the measurements [2lj. 

On the other hand, the entropy, 5", counts the (log- 
arithm of the) number of microscopic configurations or 
the multiplicity of a configuration. In the Ising model 
this is the number of configurations of spins that lead to 
the same energy. The change of entropy from a state 
with energy Ui to another with energy U2 is the number 
of ways in which spins can be flipped to reach the new 
state, i.e. the number of ways the same energy differ- 
ence, U2 — Ui, can be achieved. Let us define an anal- 
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FIG. 3: Graphical comparison of the Ising model and the one-way computer. The upper graphic, a), shows the formation 
of an ordered (magnetic) phase of spins with decreasing temperature in the 2D Ising model. The lower graphic, b), shows 
the evolution of the 2D cluster computer with increasing time to the ordered state of completely disentangled qubits. When 
this state is reached a particular solution of the computation has been chosen and can be read out from the qubits by simple 
measurements. The different stages during these processes are shortly summarised here, a) conf. 1: Ti ~ big: Spins are 
randomly flipped so that half of them point up and half of them down, without any kind of order. The magnetisation is zero. 
Ml — 0. conf. 2: T2 ^ Tcrit: Spins in the lower half and upper half point in opposite directions and form a Weiss-domain 
each with distinct magnetisation. Nevertheless, the total magnetisation is still zero. conf. 3: T3 — » 0: All spins point in 
one direction and form one ordered phase with M3 — Mmax. The critical temperature, Tcrit, for this phase transition is passed 
somewhere in between conf. 2 and conf. 3. b) conf. 1: ti = 0: Start of computation with full cluster state. All qubits 
are equally entangled and no bias towards any solution is taken, Ii — 0. conf. 2: i2 > 0: A few measurements have been 
made, some qubits are therefore not entangled with the cluster anymore, however no useful information has become available 
yet. conf. 3: ta — > tend: AH qubits have becomes disentangled and ordered and their value can be read out. The solution 
of the computation is thus fully available and the information is I3 = /max. The critical time, tcrit, at which we start to gain 
knowledge about the solution of the computational problem is passed somewhere in between conf. 2 and conf. 3. 



ogous quantity for one-way computation, the computa- 
tional capacity, C, that counts the number of configura- 
tions of qubits with the same amount of entanglement, 
E. The change of computational capacity, AC, is again 
the number of ways in which the same diiTerence in en- 
tanglement can be achieved. Due to the non-reversible 
nature of the one-way computer, a direction in time is set 
and any state with entanglement Ei will always evolve to 
a state with less entanglement, E2 < Ei. The change of 
computational capacity AC = C2 — Ci when measuring 
and transforming a state with Ei into a state with E2 is 
the number of possible evolutions (c.f. Toffoli's definition 
in [l|) that lead from Ei to £'2- The number of different 
evolutions of the computer is just the number of differ- 
ent computations, hence AC also expresses the capacity 
the initial state has to compute for a given reduction of 
entanglement. In analogy with thermodynamics and to 
allow for an additivity law for combining systems similar 
to the ones in statistical mechanics we define the com- 
putational capacity to be the logarithm of the number of 
ways in which the entanglement can be reduced. 

We note however, though the amount of contained en- 
tanglement can be calculated in symmetric multi-qubit 



states [5|, the evaluation becomes very difficult for non- 
symmetric states. Thus, it is not clear how to specifically 
calculate the computational capacity for a partially mea- 
sured cluster state. Yet the hope is, that at least bounds 
on the entanglement can be derived that are based on 
the symmetric portion of the remaining cluster. Hav- 
ing established the similarities of energy and entangle- 
ment, and entropy and computational capacity, what are 
the equivalents of the other parameters that feature in a 
phase transition? Clearly, when cooling the Ising model 
the temperature is the process parameter that is contin- 
uously reduced. How does the temperature come in and 
what is its analogue in the one-way computer? 

The set of states of a thermal system that balance the 
energy and entropy to minimise the free energy can be 
described by a curve in the set of all states, parametrised 
by a single parameter, the temperature, T. These equi- 
librium states are the thermal states. In the Ising model, 
the temperature can also be understood as setting the 
scale of how costly a perturbation from the ordered state 
is in terms of entropy. In one-way computation the state 
of the computer is changed by us as we measure the clus- 
ter to drive a certain computation forward. The pro- 
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duced curve of states that is generated by these succes- 
sive manipulations is again a curve in state space and can 
be parameterised. Indeed, since the time t moves for- 
ward monotonously during the computation we choose 
the time as the parameter of the curve or, for reasons 
that will become apparent later (c.f. Fig. [3|), wc choose 
i. The inverse time, j, tells us how much each computa- 
tional step costs. Yet, t is not simply proportional to the 
size (iV) of the cluster state since one can perform more 
than one measurement in a single computational step. 
In fact, all measurements belonging to an equal-time set 
Qt can be done in one instance of time (by definition of 
the sets Qt). Thus, t defines a measure of the depth of 
the computation and will be used to count the number of 
equal-time sets needed to successfully obtain the solution 
to a computation. 

Finally, the central role in a phase transition process is 
played by the free energy which must be balanced across 
the phases. In general thermodynamics, the free energy 
determines how much of the total energy, U = F + T S, is 
available to perform work and what portion, T S, needs 
to remain to maintain thermal fluctuations at a given 
temperature. A state of thermal equilibrium is reached 
when the free energy is minimal. To achieve the minimal 
free energy a system may take on qualitatively differ- 
ent configurations, the phases, which alter the use of the 
available energy. This qualitative change of configuration 
is what we call a phase transition and the temperature at 
which this happens is the critical temperature. For the 
one-way computer we define the computational potential 

P^E-^C, (10) 

as the quantity that balances the available entanglement 
with the computational capacity of the remaining re- 
source state at each time step in the computation, much 
like the free energy in Eq. (O balances the energy and 
the entropy for each temperature. Wc propose that the 
best way to run the one-way computer is to always min- 
imise the computational potential during computation. 
Apart from the analogy to the phase transition process 
the reason is the following. 

Universality in the context of the one-way cluster com- 
puter means that for any computational task there exists 
a cluster state of N qubits, with N big enough to pro- 
duce the solution to the computational task. Measuring 
a cluster state in such a way that the largest number 
of evolutions (or computations) remains, leaves the clus- 
ter computer most 'universal'. So, any computational 
process should be run in such a way, that the number 
of possible evolutions, the computational capacity C, of 
the remaining resource state, at each step in the com- 
putation, is maximised for a given reduction of the en- 
tanglement, E. Maximising C for a given E and time 
corresponds exactly to a minimisation of the computa- 
tional potential P. This is analogous to thermodynamics 
where the maximisation of the entropy, i.e. the maximi- 
sation of choices, for a given internal energy leads to the 



minimisation of the free energy. Of course, computations 
could be run in a less resourceful way, but for our argu- 
ment it is sufficient to discuss the optimal strategy. This 
is our proposed law of least action for the dynamics of 
the one-way computer. In analogy to the free energy, the 
computational potential thus measures how much of the 
entanglement in the resource state is 'free' to be used 
for the computation while keeping enough structure of 
connecting entanglement. 

We have identified general thermodynamic quantities, 
such as the internal energy, the entropy, the temperature 
and the free energy, that feature in all thermal processes. 
We have also established a law analogous to the second 
law of thermodynamics that makes use of the resource 
state in the optimal manner. Nevertheless, in a phase 
transition processes an additional parameter is needed 
to describe the system fully. This order-parameter ap- 
pears when a new, ordered phase emerges, for instance 
the magnetic phase in the Ising model. The parameter 
in this example is the magnetisation M and it specifies 
in which spatial direction the magnetisation points. In 
other words it serves to specify how the magnetic mate- 
rial has broken the symmetry of the governing Hamilto- 
nian. In one-way computation, 'symmetry' is broken as 
soon as we start to learn something about the solution 
to the computation. Different solutions of the computa- 
tion are like the different choices of the directions for the 
magnetisation. The order-parameter for one-way com- 
puting, analogous to the amount of magnetisation in the 
Ising model, is then proportional to the number of re- 
trieved bits of information of the solution, /. Since the 
measurements in the first stage of one-way computation 
are just for preparation purposes, they do not give us 
any information. This implies that no 'choice' has yet 
been made and the information about the solution re- 
mains zero, I{t) = for t < tcrit (early times ^ high 
temperatures). However, at some time icrit a- part of 
the former cluster state may already be deterministically 
available for readout, i.e. some qubits have become dis- 
entangled from all other qubits and are either in state 
|0) or in |1) while other qubits are still in superposi- 
tions making the results of their readout measurements 
probabilistic. Finally, when the proper read-out stage is 
reached at tend the full outcome is available and the infor- 
mation becomes maximal, /(tend) = /max- The growths 
of information over increasing time in the computation 
thus resembles the rise of magnetisation with decreasing 
temperature in the Ising model and the two curves are 
displayed schematically in Fig. [H 

Peierls' argument for the cluster computer Hav- 
ing identified the computational characteristics that cor- 
respond to the thermodynamic quantities featuring in a 
phase transition process, we now return to Peierls' argu- 
ment and apply it to one-way computation. Let us re- 
gard the final state of the computation, where all qubits 
have been measured and all entanglement has vanished as 
the ordered state of the computation. Actually, what has 
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FIG. 4: Order parameters in the 2D Ising model (magnetisation AI{T)) and the 2D cluster computer (information I{t)) are 
sketched over their process parameters, temperature T and time t, respectively (they are inversely proportional as indicated 
in the text). In 2D there exists a finite T (t) at which the order parameter M (I) becomes non-zero. In the cluster state 
computation in the first stage {t < tcrit) all measurements are performed in the equatorial x — ?/- plane to drive the computation. 
The measurement outcomes have always random outputs, i.e. the probability for either of the two outcomes is always 1/2. 
Outputs remain random until the computational phase reaches fcrit, at which the first qubits become disentangled from all 
others and their value is deterministically available. The actual outcomes are then worth one bit of information for each 
measured qubit. To obtain the full solution however, the full algorithm has to be implemented which is finished only at tend- 
The total available information grows with time to the maximal information of the outcome /max at tend- 



become ordered is the deterministic state of each single 
qubit in the cluster that implies the full knowledge of the 
solution of the computation. We can now ask the same 
question as in the case of the Ising model: Is the ordered 
state stable, or in the language of computation, will the 
solution ever be reached and the computation comes to 
an end? Again, consider first the ID case of a chain of 
qubits. The final state of the cluster computer can be 
depicted schematically as below. 

final state •••••• 

The smallest perturbation that disturbs the order is a 
single pair of entangled qubits, which leaves the last bit 
of information unavailable, i.e. the phase of completely 
disentangled, ready to be read-out qubits is interrupted. 

origins • • • • — • •, • • • — • • •, 
• • — • • • •, ... 

All states of this kind form the set of possible origins 
before reaching the final, ordered state. Again, like in 
Peierls' argument we find A^ different evolutions towards 
the ordered state. This implies an increase of the compu- 
tational capacity under the perturbation by AC = InA^. 
The change of entanglement on the other hand is just 
one unit of entanglement, AE = 1, and hence the change 
of the computational potential is AP = 1 — j\nN- The 
change should be positive for some time t to allow the 
final step of the computation to actually be executed 
at that time. However, the number of qubits in the 
cluster, A^, can be arbitrarily large to allow more and 
more complex computations. For these increasingly big- 
ger cluster states the change of the computational poten- 
tial is never negative, unless the time is infinitely large 
{t (xlnN oo). In other words, more complex compu- 
tations will not reach the final solution in a finite time 
and the one-dimensional cluster states are therefore not 
universal. 



But what is the reason that the final state cannot be 
reached? Let us take a look at the beginning of the com- 
putation to get an idea of what goes wrong. Initially we 
have the fully symmetric cluster state and we start the 
computation by measuring at least one qubit. The num- 
ber of ways in which this can be done is A^ again and the 
change of the computational capacity is InA^ while we 
have reduced the entanglement by only 2 units. Physi- 
cally this means that we have lost far too many possible 
computations (In A^) compared to the number of compu- 
tational gates that have been executed. While this im- 
balance may still permit us to solve a limited number of 
computations, it is not enough to achieve universal com- 
putation, even when we let the size of the initial cluster 
state grow larger and larger. To summarise, what is im- 
portant for one-way computers is the trade-off between 
the entanglement that they contain and the computa- 
tional capacity. The former represents the actual gates 
that our computer is able to execute, while the latter is 
the number of all possible computations that we can in 
principle perform. In ID, the change in computational 
capacity by far outweighs the number of gates imple- 
mented and this is why we cannot perform any computa- 
tion we desire. The cluster computer gets stuck at some 
point. Following the same argument one can show that in 
two dimensions universality is at least possible because 
the two terms in the computational potential grow at 
the same rate. Indeed, it is known that 2D clusters are 
universal from a straightforward constructive argument 
0. In the last section we will now quantify how the two 
quantities can be balanced in higher dimensions and with 
longer ranged correlations. 
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V. IMPLICATIONS AND CONCLUSIONS 

A number of implications are generated by the con- 
jectured analogy. Here we choose to focus on two main 
ones. First, we look at how the efficiency of computation 
scales with the dimensionality of the computer. Secondly, 
we look at the effect of allowing long range C-phases to 
generate entanglement between distant qubits while still 
remaining in ID. 

Let us calculate the change in the computational po- 
tential for a computer in d spatial dimensions. Then 
each qubit has at most 2d nearest neighbours. Note, 
that a low-dimensional system can also effectively be- 
come higher dimensional by increasing the intcrconncc- 
tivity between the qubits via non-nearest neighbour C- 
phase gates. We now repeat the argument of Peierls' for 
general d dimensions. After little inspection, it is clear 
that the optimal Peierls' cut reduces entanglement by 
an amount proportional to dN. At the same time, the 
number of possible possible cuts to achieve the same con- 
figuration is proportional to Ind^ . Balancing the two, to 
maximise the computational potential we obtain 



Ind^ cx dN, 



(11) 



giving us the critical time for computation of the form 



^crit 



Ind 



(12) 



This prediction agrees with our intuition that computa- 
tion should be faster with increased dimensionality be- 
cause there is more connectivity inside the cluster in 
higher dimensions. 

Secondly we turn to the issue of how to make a one- 
dimensional cluster computer universal. It is clear that 
we have to increase the intcrconnectivity between the 
qubits in some way. We can do that by allowing C- 
phase gates to act not only between nearest neighbours, 
but also between more distant qubits, i.e. the C-phase 
gates become long range. Therefore Peierls' argument 
has to be reviewed for the case of long range interac- 
tions. This was done in 1969 by Thoulcss in [15| and the 
result, within our context, is that if C-phases drop faster 
than the square of the distance between the qubits, then 
we cannot achieve universality of computing. A slower 
drop with distance could possibly be sufficient to allow 
universal computation. Any successful one-way model 
architecture will have to take this issue of connectivity 
into account. 

In this paper we have argued that one-way quantum 
computation can be viewed and understood as a phase 



transition. In standard phase transitions temperature is 
the parameter which determines which phase the system 
is in. In computing, this role is played by time. How- 
ever, it is transparent that temperature also has a role 
to play in computation. After all, should the ideal clus- 
ter computer be exposed to an environment of non-zero 
temperature, the thermal fluctuations would destroy the 
perfectly symmetric arrangement of exact C-phase gates. 
The higher the temperature of the environment, the more 
mixed the states of the cluster computer will be. At 
some high enough temperature this mixing will destroy 
so much of the available entanglement that the computer 
will lose all its quantum computational power. (For a 
discussion of the effect of thermal fluctuations in cluster 
computers, see 0.) 

Is it possible for us to estimate this temperature be- 
yond which the one-way quantum computer becomes 
classical? We know that in many-body systems entan- 
glement exists only below some critical temperatures. 
A very general rule of thumb is that the temperatures 
should be lower than the strength of coupling between 
the individual systems for the whole systems to be entan- 
gled [l^. An interesting conclusion can now be drawn 
when this result is translated to cluster computing. Can 
we say that the cluster state is only able to be a univer- 
sal quantum computer at temperatures smaller than the 
coupling constant for generating the C-phases? 

We leave this, along with many other issues, as an 
open question. Our conclusion is that the analogy be- 
tween thermodynamics and computation allows the flow 
of ideas between two established theories and promises 
mutual benefit. While some of the raised questions may 
be difficult to resolve, we hope that our analogy will stim- 
ulate further research into this exciting area. 
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